Effective Hamiltonians in quantum physics: resonances and geomètric phase 
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Effective Hamiltonians are often used in quantum physics, both in time dependent and time 
independent contexts. Analogies are drawn between the two usages, the discussion framed 
particularly for the geomètric phase of a time-dependent Hamiltonian and for resonances as 
stationary states of a time-independent Hamiltonian. 
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I. INTRODUCTION 

Effective Hamiltonians can arise in a variety of con- 
texts. Chosen to focus on some particular aspect or sub- 
system, an effective Hamiltonian ií c ff is constructed from 
the full Hamiltonian H of a physical system. Often, H e g 
is simpler and of dimension smaller than H . The use of 
complex energies, with the imaginary part a stand-in for 
the usually complicated and infinite-dimensional aspects 
of friction or dissipation, is a familiar examplc, occurring 
in various areas of physics 1 . Both the time-dependent 
and the time-independent Schròdinger equation admit 
dcscriptions in terms of effective Hamiltonians, forming 
the theme of this article. 



II. TIME INDEPENDENT CASE 

Consider first a time-independent situation, namcly, 
stationary states of a time-independent Hamiltonian H . 
These may include discrete bound states with negative 
energy eigenvalues and scattering states at positive en- 
ergy but also resonances in that positive energy sector, 
which are "quasi-bound" states 2 . Resonances are typi- 
cally viewed in terms of a decomposition H = Hq + V, 
and the basis states provided by the eigcnfunctions of 
Hq. These will include both bound and continuum states. 
The "interaction" V, which is contained within the full 
H, mixes these states of Ho, such a superposition con- 
stituting the quasi-bound resonance state of the full H. 
Because of this superposition, the resonance state has 
both discrete and continuum character. Further, the in- 
clusion of the continuum means that the superposition 
necessarily embraces an infinity of states. 

The splitting H = Hq + V may be either in real con- 
figuration space or in state space. The former can arise, 
for instance, when the potential in a two-particle sys- 
tem goes asymptotically to zero so that for any positive 
energy the partides separate to infinity, but intervening 
barriers may temporarily trap the system before quan- 
tum tunncling allows escape. This last phrase with its 
"before" and "temporarily" introduces a time aspect into 
a tinic-indcpcndcnt problem. While unnecessary, this ex- 



cmplifies the role of complementary pictures in quantum 
physics. The problem as a whole is of time-independent 
stationary states, some of which are characterized not 
just by an energy position but also a width (tunneling 
or otherwise) and, possibly, othcr real parameters as well 
(the so-called "profile index" q of an asymmetric reso- 
nance being an example 3 ' 4 ). The energy and width may 
be subsumed into a complex energy E r + iT/2, although 
there is nothing intrinsically complex about the prob- 
lem of stationary states of a real H . The width F in a 
time-independent picture is complementarily related to 
the lifetime of the resonance in terms of a time evolu- 
tion. Examplcs of potentials with intervening barriers 
are many: alpha-decay of nuclei, elèctric field ionization, 
two-valley potentials with angular momentum barriers 
in àtoms and molècules 3 , etc. See Fig. 1. In any of 
these systems, if the radial variable r (or a hyperspheri- 
cal equivalent in a many-particle system 3 ) is split into an 
inner and an outer region so that the former includes the 
barriers, even the states of positive E r are bound states 
within that region. But they leak out into the outer re- 
gion in the full problem and are, indeed resonances. Be- 
cause of the key role played by shapes of potentials that 
give rise to them, they are called "shape" resonances. 

Instead of the splitting in the above paragraph into two 
parts in r, the H = Hq + V breakdown may be in terms of 
states. Thus, consider doubly-excited states in àtoms and 
molècules, the first of these in the helium atom serving 
as a good example 3 . As in Fig. 2, consider states of this 
two-electron system at energies around 58 eV above the 
ground state of the atom. In 1 S e symmetry, these being 
good quantum numbers of the whole system (we will con- 
sider the non-relativistic H, ignoring spin-orbit or other 
rclativistic complications), there lie here the states of the 
one-electron ionization continuum built on the He + ls 
iònic ground state. In an independent electron descrip- 
tion, that is, when Hq includes only one-electron terms 
and possibly a mean field of the electron-electron interac- 
tion, with the residual part of this interaction constitut- 
ing V, such states can be designated lsfcs 1 S' e , where k is 
the wave-number of the continuum electron related to its 
kinetic energy through E = (hk) 2 /2m. The vàlues of E 
are, approximately, (58-24.6 = 33.4 eV). But there also lie 
in this region bound states of Hq that may be described 
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FIG. 1: Potentials with intervening barriers leading to shape 
resonances with their quasi-bound wave functions shown. (a) 
alpha-decay of nucleus of charge (Z + 2) with schematic inner 
nuclear attraction and Coulomb repulsion of separating par- 
tides, (b) Field ionization from metals in an elèctric fïeld F 
of states near the Fermi energy shown as zero. (c) Two valley 
potentials in àtoms and molècules formed from superposition 
of angular momentum and atòmic or molecular potentials, 
typically for í > 2. 



as 2s 21 S e and 2p 21 S e , which are states of an electron 
bound to He + 2s or 2p that share the same overall sym- 
metry. Thus, in terms of such an independent electron 
description, all these states are degenerate and, by virtue 
of the residual V, are superposed in the physical eigen- 
states of the full H . These are the doubly-excited "Fesh- 
bach" resonances 5 , having both discrete and continuum 
character in their description, that may be seen either 
in excitation cross-sections from the ground state (with 
58 eV of excitation energy delivered by some means) or 
in elàstic scattering of electrons from He + ls around 33.4 
eV. 

In summary of the above two paragraphs, the two al- 
ternative pictures are but just that, our pictures or de- 
scriptions. The physical system is one of eigenstates of 
the three-body Hamiltonian H in a certain energy range 
and is an integral whole. However, in our handling of 
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FIG. 2: Ground and first excited state (there are infinitely 
many) of He + and the lowest bound states of He (again there 
are infinitely many) attached to them are shown, along with 
their energies below the fully dissociated limit of two electrons 
and the helium nucleus. Cross-hatched region is the ioniza- 
tion continuum of He built on the ground state of He + , that is, 
the two-electron states lsks, E — ^k 2 being the energy of the 
continuum electron. The 2s 2 and 2p 2 states lie embedded in 
this continuum, that is, are degenerate with lsks and, there- 
fore, mix with it to give the quasi-bound resonance states. 
The resonances can be accessed either by (e + He + ) scat- 
tering or by photoexcitation from the ground state as shown 
(because of dipole selection rules, a single photon will reach 
similar states of 2s2p 1 P° symmetry whereas two-photon ab- 
sorption would be necessary to reach the 2s 2 /2p 2 1 S states). 



them, and even more in the pictures we develop for 
our understanding, we develop alternative breakdowns 
in terms of simpler sub-systems, either in the real con- 
figuration space of r or in two-electron states of an inde- 
pendent particle picture. Even as we do so, pictures and 
models being necessarily involved in our approach to any 
physical situation, we must also keep in mind that Ho, 
V, state configurations 2s 2 , 2p 2 and lsks, etc, are not 
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elements of the underlying physical reality nor are they 
accessible to measurement. In the same vein, the distinc- 
tion between the terms "shape" and "Feshbach" is also 
somewhat arbitrary but, nevertheless, useful. In partic- 
ular, with reference to Figs. 1 and 2, whereas the formcr 
often lie just above the energy threshold to which they 
are attached, the latter lie just below (the He + 2s thresh- 
old). Correspondingly, the shape resonances tend to be 
broad (shorter lived) whereas the Feshbach resonances 
are typically narrow (longer lived). In a closely-related 
two-electron system to our discussion above, namely the 
negative ion of the hydrogen atom, and in the symmctry 
1 P° which can be reached by one-photon absorption from 
the ground state (unlike 1 S e which takes two photons), 
an example of each was clearly shown on either side of 
the H(n = 2) threshold by a clàssic experiment 6 . See 
Fig. 3. 

The partitioning of a Hamiltonian into two spaces is 
conveniently done in the "Feshbach projection operator" 
formalism 8 . With P and Q projection operators into 
what are usually referred to as "closed" and "open" sub- 
spaces (for alternative terminologies, see Section 8.6.1 
and 8.6.2 of Ref. 3), the former a finite space with 
bound eigenstates and the latter the infinite space of the 
continuum into which decay takes place, we can rewrite 
Hip = Eijj as 



H(P + Q)ip = E(P + Q)r/>, P + Q=l. 



(1) 



Multiplying from the left by Q and re-arranging with the 
use of Q 2 = Q, a projection property, we have 
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Qip = (E - QHQY 1 {QHP)ip. (2) 

Multiplying Eq. (1) by P from the left, and substituting 
for Qtp from Eq. (2) gives 

[PHP + PHQ(E - QHQ)- 1 QHP)]Pi/j = EP^. (3) 

The term in square brackets is the effective Hamiltonian 
H fi in P-space. The term "optical potential" is also 
used. Note that it has P at both right and left extremes. 
With tp also occurring only in its projected piece, the 
Eq. (3) is cntirely in P-space. With no approximation 
made, the second term in H e g includes the effect of the 
remaining Q-space. Note its structure in the form of 
a second-order energy, with H "carrying" from P to Q 
space and, then with an attached energy denominator, 
entirely in that Q-space, the H "returning" again to P- 
space. The first PHP term has a purely discrete spec- 
trum but some of thosc eigenvalues may acquire both a 
shift and a width as a result of the second term, which 
involves coupling to the Q-space. These are the reso- 
nances. Because of the energy denominator, the effect 
of this second term can be dramàtic at energy vàlues E 
close to an eigenvalue of QHQ. 



FIG. 3: Theoretical hyperspherical potential wells and exper- 
imental observations of doubly-excited P° resonances of H~ 
in the vicinity of the H(iV = 2) threshold as observed in pho- 
todetachment. The cross-section (in units of squared Bohr 
radius) shows a sharp Feshbach resonance just below and 
a broad shape resonance just above the 10.95 eV threshold 
(from Ref.6). The potential wells (1 Ry =13.6 eV) in hyper- 
spherical coordinates for the two-electron system account for 
the observed resonances, the barrier in the + curve holding 
the shape resonance while the — well supports the Feshbach 
resonance 7 . 

III. TIME DEPENDENT CASE 

Turn next to the time-dependent Schròdinger equa- 
tion, iU(t) = H(t)U(t), with a dot standing for a 
time derivative, and where we choose to work with 
the evolution operator U rather than the wave func- 
tion ip(t), to which it is easily related through ip(t) — 
J7(í)V>(0), 17(0) = 1. Upon writing 

U(t) = U 1 (t)U 2 (t), (4) 

taking a derivative with respect to time, and multiplying 
from the left by Uy , we get 
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iÜ 2 = H eff U 2 , H eff = U^HUi - íU^Ül (5) 



M2 = M y(l + | M3 | 2 ), e 1 ™^ =(1 + | M3 | 2 ), (9) 



Again, we have a reduced equation for U 2 alone but, of 
course, thc cfïcctivc Hamiltonian incorporates the part 
containcd in U\ so that the reduced expression is formally 
complete with no approximation implicd. A connection 
to variational principies and identities will be discussed 
at the end of this Comment. 

The above is completely general for any time- 
dependent problcm. One important specific applica- 
tion is in separating the geomètric phase 9 from other 
elements of time evolution. With the advent of fault- 
tolcrant quantum computation, the geomètric phase is 
seen to have advantages in robustness and fidelity over 
dynamical phases that accumulate as a result of energy 
changes in time 10 . Consider, as an illustration, an angu- 
lar momentum that couples through its magnètic mo- 
ment to a magnètic field B(t), with only linear cou- 
pling in H = —J-B(t). Given the three operators 
in the system, conveniently chosen as the usual triad 
( J + = J x + iJ y , J_ = J x — iJ y ,J z ), a complete solution 
for the evolution operator can be obtained by writing 11 



U(t) = e -iMt)J+ e -iMt)J- e - 



im(t)J z 



(6) 



That this is indeed the solution can be seen by a con- 
struction that leads to the required equations that define 
the functions // in the exponents. Taking the time deriva- 
tive, rcpeatcdly applying a Standard, "Baker-Campbell- 
Hausdorff" (BCH), identity 12 for e A Be~ A to cast as an 
operator multiplying U from thc left, gives 



ÍÜ = [fl 3 J+ + (Jz + Íp 3 J+)fll 

+ (/i 2 - iH2Íi\){J- - ZilJ-sJz + H3 2 J+)]U- (7) 

Upon identifying the right-hand side of Eq. (7) with HU, 
the equations satisfied by the three \x follow: 



M3 + ^M3 2 -B+ - i/j, 3 B 3 = ~B_, 

fi 2 - = ~\ B +> í" 1 + iB +^ = ~ B s, (8) 

where B± = B x ± íB y . The requirement U(0) = X sets 
the boundary conditions for all three functions, /i(0) = 0. 

Thus, solutions of the set of classical equations in 
Eq. (8), when inserted into Eq. (6), give the full quan- 
tum evolution. Thc first equation for ji + decouples from 
the other two and may be solved by itself first, followed 
by simple quadrature of the other two. Note that any 
value of spin- j obeys the same set of equations in Eq. (8), 
since no use was madc of any specific represcntation of 
the operators J, only their commutators. Further, U is 
guaranteed to be unitary by construction 11 . Indeed, this 
implies relationships between the three complex fi: 



so that there are only three linearly independent quanti- 
ties which may be chosen as the real and imaginary parts 
of yu 3 and Re /Ui. With this, for j = 1/2, Eq. (6) may be 
written as 



e-^/ 2 ( 1 -i/ise^" 1 



U(t) 



Thus, an initial state ip(Q) = ( q I cvolves to 



(10) 



m = 



-ÍMl/2 



1 + M 2 V -*/*3 

and the density matrix p(t) = \ip(t)){vp(t)\ becomes 



(11) 



P{t) 



1 



|A*3| 



1 



-*M 3 lM3| 



(12) 



Whereas the density matrix involves only the two param- 
eters of the complex quantity yti 3 , and so does the wave 
function (except for a usually unobservable phase /xi), 
the evolution operator depends on the third parameter, 
Re pi, as well. The above equations give the 2x2 matri- 
ces for a spin-1/2 but all these features of the role of thc 
three parameters contained in /i apply also to thc (2j + l) 
vectors and matrices of any spin- j. 

The phase pi in the above expressions, particularly in 
the evolution operator U in Eq. (6) may be viewed as 
an illustration of Eq. (4), where U\ is the product of the 
first two exponential terms in Eq. (6) and Ui the last 
term involving \i\ alone. We have, upon using Eq. (9), 



fi 



zv4 2 



(13) 



which depends only on fi 3 whereas fi\ is contained in 



Un = 



e -Vl/2 







oVl/2 



(14) 



It is instructive to see the specific form, including of in- 
dividual terms, that H e g in Eq. (5) takes in this examplc. 
We have U^HUi = 



i l· i* 3 B--i l· L3B + -B 3 (l-\n3\ 2 ) ■ r> 1 R l.,2 R 

2(i + | M3 |2) ^3^3 - 2 tí - - 2^ tí + 

-2i l·í lB 3 -B + -\ f i 3 \ 2 B- i l·í3 B + -i f i*B-+(l-\ l·í3 \ 2 )B 3 



and iU 1 1 Úi = 



(15) 



■5 



i^B_+^3|M3| 2 -B++2|M3| 2 B 3 

2(i+ImI 2 ) 

-2i At ;5 3 -B+-|M3| 2 -B- 

2(TT^F? 



- t il(i t ilB++iB-+2 l· L 3 B 3 ) 



(16) 

Upon subtracting Eq. (16) from Eq. (15) to form H e g 
in Eq. (5), the off-diagonal terms cancel. The diagonal 
terms add to give precisely (— B3 — i^i 3 B + )a z /2 to coin- 
cide with the last Eq. (8) which gives the phasc fj,i. Our 
analysis for the triad choice (J+, J_, J3) in Eq. (6) could 
also have been carried out for other choices of three lin- 
early independent operators, for example, the Cartesian 
(J x , J yi J z ) or an Euler-angle set. 

Having provided the spin-1/2 or SU(2) group's decom- 
position of the cvolution opcrator into two factors U\ 
and U2 in detail as a pedagogical illustration, wc note 
that we have given a closely paral·lel development for 
a more general SU(iV) for an arbitrary time-dependent 
Hamiltonian 13 . This construction reduces inductively 
the operator for N to the one for (N — 1) with dcfining 
equations for complex parameters z(t) which are analogs 
of (U 3 and phases that are analogs of Further, while 
the U (t) in Eq. (6) is always unitary by construction, the 
individual U± and U2 above are not but can also be unita- 
rized. This is accomplishcd above by separating fií in U2 
in Eq. (14) into real and imaginary parts, the latter as in 
Eq. (9) incorporated into U\ to make it unitary, leaving 
U2, which depends on Re Hi, as a pure phase. A very 
similar construction 13 for general SU(iV) also provides an 
explicitly unitary decomposition of the evolution opera- 
tor, with (N — 1) phases, and again their decomposition 
into dynamical and geomètric pieces as in Eq. (15) and 
Eq. (16). 



IV. RELATION TO VARIATIONAL 
PRINCIPLES AND IDENTITIES 

Finally, the connection of the effective Hamiltonian 
in Eq. (5) to variational principies and their associ- 



ated identities 14 is also instructive. Following a general 
construction , the time-dependent equation ill = HU 
can be converted into an identity 



U(t) = U t (t) - / dt'L(t')[iU t (t') - H(t')U t (t% (17) 
Jo 

where U t is a "trial" function and L a "Lagrange adjoint" 
function [14] given by 



L(t') = -iU(t)U-\t'). 



(18) 



The identity in Eq. (17) is easily verified upon doing an 
integration by parts and using Eq. (18). Combining the 
two equations, we have 



U(t) = U t (t) + iU{t) / dt'U~\t')[iÜ t (t') - H(t')U t (t% 
Jo 

(19) 

and an associated variational principle for U val which fol- 
lows upon replacement of U on the right-hand side of 
Eq. (19) by its trial approximation U t - 



U yBI {t) = U t (t)(l-i í dt'Ut\t')[-iÜ t {t')+H{t')U t {t')]). 
Jo 

(20) 

The right-hand side can be evaluated with any approxi- 
mate solution Ut and, as is clear from the derivation, will 
give a variationally correct approximation with only sec- 
ond order errors in (U t — U). Thus, U val improves on U t 
which is only correct to first-order. Note the appearanec 
of the H c s of Eq. (5) in the integrands of Eq. (19) and 
Eq. (20). 
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